Energy transport in a one-dimensional granular gas 
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We study heat conduction in one- dimensional granular gases. In particular, we consider two 
mechanisms of viscous dissipation during inter-grain collisions. In one, the dissipative force is pro- 
portional to the grain's velocity and dissipates not only energy but also momentum. In the other, 
the dissipative force is proportional to the relative velocity of the grains and therefore conserves 
momentum even while dissipating energy. This allows us to explore the role of momentum conser- 
vation in the heat conduction properties of this one-dimensional nonlinear system. We find normal 
thermal conduction whether or not momentum is conserved. 
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I. INTRODUCTION 

Energy transport in low-dimensional systems can be 
pathological in the sense that Fourier's law for heat 
conduction might break down [l|, d, H, In three- 
dimensional solids, the energy (heat) transported is gov- 
erned by Fourier's law, which says that the heat flux J 
is proportional to the gradient of the temperature, 

J = -kVT, (1) 

where k is the thermal conductivity. The thermal con- 
ductivity in Fourier's law is independent of system size 
and of time. This equation assumes that a local equi- 
librium is established at each time, so that one can de- 
fine the local energy flux J(x, t) and temperature T{x, t), 
when the temperature gradient lies along the x direc- 
tion (3]. In one-dimensional systems one might expect 
a similar equation, with the gradient replaced by the 
derivative of T(x, t) with respect to x, and J simply 
being the energy flux per unit time. However, it has 
been observed that in many one-dimensional models the 
thermal conductivity varies with the system size N as 
K - A^", with a lying between 0.32 and 0.44 . Self- 

consistent mode coupling analysis for one-dimensional 
nonlinear chains 0] predicts a universal value of a ~ 2/5, 
while rcnormalization group analysis based on hydrody- 
namical models supports a universality class with a ~ 
1/3. Recently, Mai et al. [1, @| revisited the problem and 
obtained the same a = 1/3 behavior for one-dimensional 
nonlinear chains as obtained for the fluid- like systems. 
They argue that the dynamic equations are those of a 
fluid at all length scales even if the static order of the 
chain extends to very large N. They explain the dis- 
agreement with most numerical results as a consequence 
of the numerical difficulties in reaching asymptotic be- 
havior, which requires extremely large N, but their con- 
clusions have in turn been disputed Whether there 
is a single universality class or there are two, there is a 
connection between heat conduction and diffusion in one- 
dimensional systems 0, that is, normal diffusion leads to 
normal heat conductivity (k ^ N'^)-, whereas anomalous 



transport is associated with k ^ N", with a > (< 0) 
for super(sub)diffusion. Thus, superdiffusion is associ- 
ated with a divergent thermal conductivity. 

Another key concept in the understanding of the ther- 
mal conductivity in Id systems is momentum conserva- 
tion [3]. In fact, momentum conservation usually irnplies 
the divergence of the thermal conductivity in Id [l|, 0], 
and yet it is not necessary for the occurrence of anoma- 
lous diffusion [To| . An exception to the momentum con- 
servation rule is a chain of particles interacting via the 
nearest-neighbor potential V{qi — Qi-i) = 1 — cos{qi — 
Qi-i), the so-called rotator model [11|. This momentum- 
conserving model exhibits normal transport behavior to 
very high numerical accuracy, and is said to occur be- 
cause the rotator model "cannot support a nonvanishing 
pressure, and thus infinite-wavelength phonons cannot 
carry energy" [l3| . In ^ this behavior is ascribed to the 
periodicity of the potential and the associated indepen- 
dent jumps from valley to valley. 

In this contribution we discuss heat conduction in one- 
dimensional granular gases. Granular gases are dissi- 
pative systems, the associated energy loss usually being 
modeled by introducing a coefficient of restitution as a 
parameter in the description of granular collisions. In- 
stead, we follow a more dynamical approach and intro- 
duce energy dissipation explicitly via viscous terms in the 
equations of motion for the granules. There are a number 
of different sources and descriptions of viscous effects in 
the literature [H, M, Q M, [H, [13, The different 
ways in which they appear in the dynamical descriptions 
suits our particular focus of interest, which is the role 
of momentum conservation in the heat transport pro- 
cess. In particular , on e way to introduce viscosity con- 
serves momentum [1^ [g^, [111, [13] , while the other does 
not [H, [13, [13| . Interestingly, we establish that thermal 
conduction exhibits normal behavior in both cases, that 
is, that there is no divergence as the system size increases. 

In Sec. |TT]we describe the two dissipative models, and 
in Sec. IIIII we present simulation results to characterize 
the transport of heat through the granular gas in both 
cases. We conclude with a short summary in Sec lIVI 
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II. DISSIPATIVE MODELS 

We consider N identical granules constrained to move 
on a line between two walls at different temperatures. 
The granules move freely except during collisions either 
with the walls or with one another. The system is a one- 
dimensional "granular gas" because the distance between 
the walls is much greater than the space occupied by the 
granules. The walls act as heat baths, that is, whenever a 
granule collides with a wall at temperature T, its energy 
is absorbed by the wall and it acquires a new velocity 
away from the wall according to the probability distribu- 
tion il 



Piv) = y exp {-v'llkBT) 



(2) 



Intcrparticle collisions are governed by the power-law po- 
tential 



Vi^kM^l) =0, 



(3) 



(S > 0. 



Here (S^./t+i = yk+\ — Vk, a is a prefactor determined 
by Young's modulus and Poisson's ratio, and the prin- 
cipal radius of curvature R of the surfaces at the point 
of contact [1^, [l^l ; and yu is the displacement of gran- 
ule k from its position at the beginning of the collision. 
The exponent n is 5/2 for spheres (Hertz potential), it 
is 2 for cylinders, and in general depends on geometry. 
We stress that the one-sided (only repulsion) granular 
potential even with n = 2 is entirely different from a 
two-sided (repulsion and attraction) harmonic potential. 
The one-sidedness of the potential leads to analytic com- 
plexities even in the dissipationless case [U, [2^ [29| . 
and even greater complexities in the presence of dissi- 
pation [H, 111, HO, HI . 

In this paper we explore the low density limit, which 
leads to enormous analytical and computational simplifi- 
cations. The low density feature of these approximations 
is implemented via the assumption that the collisions are 
always binary [sot , that is, that only two granules at a 
time are members of any collision event, and that at any 
moment of time there is at most one collision. 

Our approach starts with the equations of motion of 
the particles labeled by index fc, and so we write yk as a 
function of time t. It is convenient to deal with scaled 
position and time variables Xk and i, related to the un- 
sealed variables yk and r as follows, 
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Here m is the mass of the granules, and the velocity vq is 
an arbitrary choice in terms of which other velocities are 
expressed. We also introduce a scaled friction coefficient 



7 
mvQ 



l/r, 



(5) 



where 7 is the friction coefficient for the unsealed vari- 
ables. In the low density limit we only need to con- 
sider the equations of motion for two colliding granules 
A: = 1,2. Furthermore, in this paper we only consider 
cylindrical grains, which leads to considerable simplifica- 
tion while still capturing the important general features 
of the system that we seek to highlight. We stress that 
the one-sided granular potential (i.e., one with only re- 
pulsive interactions) even with n = 2 is entirely different 
from a two-sided harmonic potential. 

Consider a viscous force that is proportional to the rel- 
ative velocity of the colliding granules. Such a force has 

Jll, ill but it 



been considered not only theoretically ._ . , 

has also been observed experimentally |2ll.l22j|. While one 
might be tempted to think of this force as arising because 
one grain rubs against the other, the actual mechanism 
is more complicated and involves the medium that sur- 
rounds the granules [U, The exact mechanism is 
still a matter of conjecture. In any case, the appropriate 
equations of motion in this case are 

xi = [-j{xi - X2) - {xi - X2)]0{xi - X2), 

X2 = [7(^1 - ^2) + {xi - X2)] 0{xi ~ X2), (6) 

where a dot denotes a derivative with respect to t. The 
Heaviside function is defined as 0{x) = 1 for x > 0, 
e{x) = for X < 0, and 6'(0) = 1/2. It ensures that 
the two particles interact only when in contact, that is, 
only when the particles are loaded. The post-coUisional 
velocities (called u below) can be written in terms of the 
velocities of the two granules at the beginning of the col- 
lision (called v) as 

«i = i[l~e-^*°]z;i + i[l + e-^*°]«2. 



■e-^*°] V2, 



(7) 



where to = 7r/(-\/2 — 7^) is the duration of the collision. 
Since ui + 1*2 = wi + W2, the momentum of the center of 
mass is conserved. 

A model in which the center of mass velocity is not 
conserved is one governed by the equations of motion [l^, 

m 

Xi = [--fxi - {xi - X2)]9{xi - Xi), 

X2 = [-^X2 + {xi ~ X2)]9{xi ~ Xl). (8) 

Here one might think of the viscosity arising from an in- 
teraction with the medium. However, such an interaction 
would produce a damping term not only during a colli- 
sion but also while the granules are moving independently 
between collisions. In the low density limit the granules 
would hardly collide before stopping entirely unless the 
viscosity is extremely low, and one can then not talk of 
heat transport along the one-dimensional system. Thus, 
this model, while perhaps not realistic in any sense, is a 
"toy model" in which momentum is not conserved while 
still capable of supporting energy transport and therefore 
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relevant to our question. Indeed, we can again calculate 
the post-collision velocities, 

2 L J 2 



III. NUMERICAL SIMULATIONS 



Ul 
U2 



W2, 



(9) 

where = 27r/(y8~72) is again the duration of the 
collision. It is easy to verify that in this case the momen- 
tum of the center of mass decays from vi + V2 before a 
collision to 

Ml +U2 = bi +t^2]e-^*°, (10) 

after the collision, so that the momentum is not con- 
served. 

Caution must be exercised in the choice of the damping 
coefficient 7. If it is too large, energy acquired from cither 
wall is simply dissipated before it crosses the system, in 
which case the problem changes from one of energy trans- 
fer from one wall to the other through the granular gas 
to that of two walls pumping energy into the granular 
"sink." To estimate the limit on the damping coefficient 
we can imagine a sequence of collisions whereby a granule 
starting from the hot wall with kinetic energy collides 
with the next granule, which in turn collides with the 
next one. and so on, until the last granule, whose en- 
ergy is i?jv, collides with the cold wall. The change in 
the kinetic energy due to each collision in the momentum 
conserving model is 

^ \{ul+ul)-\{vl+vl) = -\{l- er'^'^) (-l+«2)^ 

(11) 

and in the momentum non-conserving model it is 

AE ^ ^iul+uly\ivl+vi) = (1 - e-^*«) (^^1+^.2)^ 

(12) 

In both cases the kinetic energy before and after the col- 
lision are related by an expression of the form Eafter = 
£'be/ore(l " b^tg) where 6 is a velocity-dependent dimen- 
sionless quantity of order unity, and where we have as- 
sumed that jto <C 1. As an order of magnitude estimate 
it is sufficient to write E]s[ ^ Eo{l — 7*0)^- The average 
energy of a hot granule is ksTi/i and in order for energy 
to flow across the system the energy of the granule that 
arrives at the cold wall must be greater than kBT2/2. 
Consequently, we estimate that a requirement for energy 
transport is that 

l/N 

7<l-(^) • (13) 



T2 
Ti 

We now wish to explore the following specific ques- 
tions. Is there a finite thermal conductivity for either 
or both of these models? If so, which model has a 
higher/lower thermal conductivity? How does the ther- 
mal conductivity depend on temperature? On viscosity? 
On system size? These are the questions we address in 
the next section. 
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FIG. 1: (Color online) From top to bottom: energy injection 
by the hot wall (red), energy absorption by the cold wall (ma- 
genta), energy injection by the cold wall (blue), and energy 
absorption by the hot wall (green), for a momentum non- 
conserving system of A*' = 500 particles between walls main- 
tained at Ti = 6 and = 3 and with 7 = 0.0002. The dashed 
lines are the linear regressions of these curves (after the tran- 
sient). They yield the rate of energy injection/absorption. 
The initial temperature of the granular gas is T = 1. 



The low density limit allows us to use an event driven 
algorithm in our simulations. As indicated earlier, when 
a granule collides with a wall at temperature T its en- 
ergy is absorbed by the wall and it acquires a new en- 
ergy as determined by the velocity distribution given in 
Eq. ^ . The initial temperature of the granular gas is ar- 
bitrarily set to T = 1, and we allow the system to arrive 
at a steady state before beginning our "measurements." 
While ideally one would want the temperature gradient 
in the steady state to be strictly linear as assumed in 
the linear response theory that leads to the Fourier law, 
the subtleties encountered in Id systems are well-known 
and observed pretty much no matter how the thcrmal- 
ization is implemented @, [2^ and the resulting gradients 
are only strictly linear away from the boundaries (non- 
linear behavior typically sets in near the boundaries). In 
any case, we can address the question of whether or not 
there are system size divergences in the transport of heat. 
That we do in fact achieve a steady state can be seen, for 
example, in Fig. [1] where we show the energy absorbed 
and injected by each wall as a function of time for the 
dissipating momentum model. After a short transient, 
the rate of energy injection and absorption become con- 
stant, as they should in a steady state. We ran tests 
to ascertain that the initial temperature of the granu- 
lar gas is not important for this equilibration, that is, 
we find that the rates of energy injection and absorption 
arc independent of the initial velocity distribution of the 
granular gas. The momentum-conserving model equili- 
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brates equally well, also independently of initial condi- 
tion. Henceforth we set the initial temperature of the 
granular gas to unity. 
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FIG. 2: (Color online) Rates of energy transmission as a func- 
tion of the temperature gradient. In this simulation the vis- 
cosity coefficient was set to 7 = 0.0002 and the system is 
composed of 500 granules. The (red) plus signs correspond to 
the momentum conserving system, while the (blue) stars cor- 
respond to the momentum dissipating system. In this figure, 
the cold wall temperature was set to 3. 
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FIG. 3: (Color online) Thermal conductivity as a function of 
the temperature. In this simulation the viscosity coefficient 
was set to 7 = 0.0002 and the system is composed of 500 
granules. The (red) plus signs correspond to the momentum 
conserving system, the (blue) stars to the momentum dissi- 
pating system. 

Next, we discuss the dependence of the energy flux on 
the temperature gradient. The energy flux J = dE/dt 
was calculated as the slope of the transmitted energy 
as a function of time, that is [(energy injection by the 
hot wall - energy absorption by the hot wall) - (energy 
injection by the cold wall - energy absorption by the cold 
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FIG. 4: (Color online) Thermal conductivity as a function of 
the viscosity coefficient. In this simulation the temperature 
of the cold wall is set to 3 and that of the hot wall to 6. 
The system is composed of 500 granules. The (red) plus signs 
correspond to the momentum conserving system, the (blue) 
stars to the momentum dissipating system. 
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FIG. 5: (Color online) Rates of energy transmission as a func- 
tion of the system size. In this simulation the viscosity co- 
efficient was set to 7 = 0.0001 and the temperatures of the 
walls were set to 3 and 6. The (red) plus signs correspond 
to the momentum conserving system, the (blue) stars to the 
momentum dissipating system. 



wall) - (energy dissipated during flow along the chain)] 
per unit time. As seen in Fig. [21 both models lead to 
behavior fairly well described by Fourier's law. Repeating 
this plot for different temperatures of the cold wall, we 
obtain the dependence of the thermal conductivity on the 
temperature. As observed in Fig. [31 k is an increasing 
function of T for both models, but it is larger for the 
momentum dissipating system. On the other hand, for 
a given temperature, the thermal conductivity decreases 
(almost linearly) with the viscosity (see Fig. [J). Once 
again, the dissipating system presents larger values of k. 
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FIG. 6: (Color online) Rates of energy transmission as a func- 
tion of scaled system size. The temperatures of the walls were 
set to 3 and 6. The (red) plus signs correspond to the mo- 
mentum conserving system, the (blue) stars to the momentum 
dissipating system. 



Finally, we have analyzed the dependence of the rate 
of energy transmission on system size, for a fixed temper- 
ature gradient and fixed viscosity. This is the crucial test 
of normal vs anomalous behavior. In Fig. Owe can clearly 
see that as the size of the system increases, the rate of 
energy transmission does not increase with system size. 
In fact, it decreases. Therefore, in our model the ther- 
mal conductivity does not diverge with increasing system 
size regardless of momentum conservation or dissipation. 
One may be tempted to think of this behavior entirely 
as a consequence of the energy dissipation because for 
bigger systems the number of collisions necessary for the 
energy to be transmitted from one end of the system 
to the other increases. Hence more energy is dissipated 
and less energy arrives at the cold wall. However, the 
introduction of dissipation changes the dynamics more 
profoundly. This can be seen in Fig. [51 where we plot the 
rates of energy transmission against jN, thus taking into 
account the "simple" effect of energy dissipation [see dis- 
cussion preceding Eq. p^ ] . The rate of energy transfer in 
the momentum conserving system is essentially indepen- 
dent of system size in this scaled representation, but that 
of the momentum dissipating system actually decreases 



even when scaled in this way. In any case, we note the 
very small scale of variation of the ordinate in Fig. [HI 



IV. CONCLUSIONS 

Heat transport in one-dimensional discrete systems 
continues to be a problem of theoretical interest, un- 
certainty, and even controversy [l|, H, H, 0, H @|- Very 
recently, experimental results in this arena have also 
started to appear. In (3]| the focus is the understanding 
of transport of heat along colliding granular beads in a 
liquid medium where the questions of interest involve the 
nature of the contact regions ("liquid bridges") between 
granules. In [sl] the issue is the breakdown of Fourier's 
law in nanotube thermal conductors in that the thermal 
conductivity diverges with length of the nanotube. The 
point is that even after many years of study the condi- 
tions that lead to the validity or violation of Fourier's law 
are not yet clear. While momentum conservation has of- 
ten been featured as a condition closely associated with 
the system size divergence of the Id thermal conductiv- 
ity, we have examined two dissipative granular gases, one 
of which involves momentum dissipation along with en- 
ergy dissipation, whereas in the other momentum is con- 
served. In both of these the thermal conductivity remains 
finite and in fact decreases as the number of granules in 
the system increases, as it must in a system where each 
collision leads to energy dissipation. In the momentum 
dissipating model, it decreases more rapidly than can be 
accounted for by the simply heat loss mechanism of the 
collisions, indicating a more profound change in the dy- 
namics. This behavior points to the caution that must 
be exercised when associating momentum conservation 
with anomalous behavior in one dimension. 
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